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SUMMARY 


The  purpose  o£  this  research  was  to  determine  a  rational  basis 
for  interpretation  of  the  aerodynamic  force  characteristics  of  helicopter 
rotors.  The  approach  consisted  of  three  elements.  The  first  was  the 
pa  iliel  development  of  series  expressions  for:  flight  measurements  of 
rotor  differential  pressure  distributions;  corresponding  two-dimensional 
wind-tunnel  pressure  distributions  for  the  same  airfoil  geometry;  and 
theoretical  pressure  distributions  for  the  two-dimensional  airfoil.  Second, 
there  was  a  comparison  between  pressure  differential  distribution  compo¬ 
nents  determined  from  flight  measurements  and  pressure  distribution 
components  computed  by  an  approximate,  unsteady,  three-dimensional 
theory  developed  previously  at  CAL  for  the  U.  S,  Army.  Third,  the  sen¬ 
sitivity  of  the  airfoil  drag  characteristics  to  variations  in  the  pressure 
differential  distributions  was  investigated  on  the  basis  of  two-dimensional 
steady  boundary  layer  theory. 

The  fundamental  technique  employed  was  a  curve -fitting  process 
based  on  an  expansion  of  the  pressure  distributions  in  terms  of  Glauert 
coefficients.  Available  flight  measurements  of  pressure  distributions  for 
the  UH-1A  rotor  at  advance  ratios  of  y.  =  0,  =  0.08,  and  jli-  =  0.26, 

and  the  H-34  rotor  at  advance  ratios  of  jj,  -  0.  18  and  yu,  =  0.  29  were 
used  in  the  analyses.  Corresponding  wind-tunnel  measurements  on  two- 
dimensional  NACA  0015  and  0012  airfoils  were  also  reduced  to  Glauert 
form.  There  is  only  one  coefficient  in  the  theoretical  representation  of 
the  differential  pressures  because  these  two-dimensional  airfoils  are  sym¬ 
metrical.  Induced  camber  effects  arising  in  the  available  finite  span 
theory  were  much  smaller  than  the  values  deduced  from  measured  blade 
pressure  differentials. 

Glauert  coefficients  derived  from  the  measured  rotor  pressure 
differentials  had  many  of  the  characteristics  of  the  Glauert  coefficients 
derived  from  measured  two-dimensional  pressure  distributions.  They 
also  differed,  partly  because  of  the  lack  of  a  good  quantitative  measure  of 
angle  of  attack  (which  greatly  influences  the  viscous  effects).  It  was 
found  that  the  viscous  effects  evident  in  the  two-dimensional  data  could 
also  be  seen  in  the  rotor  data.  It  is  postulated  that  these  viscous  effects 
probably  account  for  a  large  part  of  the  differences  observed  between  the 
theoretical  predictions  and  experimental  measurements  of  the  rotor 
pressure  differential  components. 

Calculated  drag  coefficient  distributions  (based  on  differential 
pressure  measurements)  exhibited  a  characteristic  rise  on  the  retreating 
side  of  the  rotor  disk. 
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INTRODUCTION 


Expansion  of  the  flight  envelope  of  the  helicopter  (and  of  air¬ 
craft  that  are  helicopter  derivatives)  is  proceeding  at  a  rapid  rate. 

Nearly  all  the  performance  gain,  however,  has  come  about  through  the 
application  of  additional  power  or  through  configurational  changes  that 
minimize  the  role  of  the  rotor  in  the  high-speed  flight  range.  Dynamic 
problems  (fatigue  of  rotating  components,  shaking  forces  transmitted  to 
the  fuselage,  etc.)  remain.  Concurrent  with  this  performance  growth 
have  been  research  efforts  devoted  to  aerodynamic  models  which  retain 
some  part  of  the  rotor  wake  effects  and  blade  dynamic  motion  effects  on 
air  loadings.  However,  both  performance  predictions  and  blade  dynamic 
response  predictions  still  are  based  largely  on  some  empirical  represen¬ 
tation  of  local  airfoil  characteristics;  in  fact,  measured  two-dimensional 
section  characteristics  are  often  assigned.  It  is  generally  believed  that 
the  use  of  two-dimensional  airfoil  characteristics  can  introduce  sizable 
errors  in  the  estimation  of  rotor  performance- -particularly  in  the  esti¬ 
mate  of  the  power  required.  Consequently,  empirical  correction  factors 
derived  from  performance  tests  are  incorporated. 

In  recognition  of  these  deficiencies  of  knowledge,  the  U.  S. 

Army  Aviation  Materiel  Laboratories  (USAAVLABS)  initiated  a  coordinated 
experimental  and  theoretical  research  program  on  helicopter  blade  loads 
(and  stresses)  several  years  ago.  In-flight  measurement  of  chordwise 
time-varying  blade  pressure  distributions  (using  the  H-34  and  UH-1A) 
was  the  major  experimental  effort;  Cornell  Aeronautical  Laboratory,  Inc. , 
(CAL)  participated  in  the  overall  program  as  a  contractor  developing 
theoretical  prediction  techniques.  A  major  advance  was  made  in  that 
accurate  data  were  accumulated,  and  good  agreement  was  obtained 
between  measured  and  theoretical  time  histories  of  section  lift  (References 
1  and  2). 


Initially,  emphasis  in  the  measurement  program  was  placed  on 
the  rotor-blade  lift  distribution;  pressure  tap  locations  and  data  reduction 
methods  were  specifically  tailored  to  minimize  the  errors  in  the  lift  deter¬ 
mination.  No  attempt  was  made  to  obtain  the  aerodynamic  pitching 
moments  from  the  experimental  data  and,  of  course,  no  direct  indication 
of  section  drag  was  obtained.  Extraction  of  additional  information  from 
the  available  measured  differential  pressure  distributions  was  the  intent 
of  the  effort  reported  herein. 

The  specific  purposes  of  the  present  research  are  (1)  to  test  in 
detail  the  current  assumption  that  the  rotor-blade  sections  experience 
pressure  distributions  similar  to  those  measured  on  an  isolated  two- 
dimensional  airfoil,  with  the  same  cross  section,  in  a  steady  wind  stream; 
(2)  to  compare  the  components  of  the  measured  pressure  distribution  with 
those  predicted  by  the  most  refined  rotor-blade  aerodynamic  theory  which 
is  readily  available;  and  (3)  to  estimate  the  profile  drag  of  rotor-blade 
elements  from  the  measured  differential  pressures  on  rotor-blade  sections. 


1 


The  measurements  are  used  in  conjunction  with  available  theoretical 
techniques  to  achieve  the  goals  listed  above.  The  central  technique  is 
the  decomposition  of  the  pressure  distribution  into  components  so  that 
quantitative  comparisons  can  be  made,  for  example,  between  flight 
measurements  on  blades  and  two-dimensional  wind-tunnel  data.  Pressure 
measurements  also  were  used  to  estimate  local  velocity  components  and 
these,  in  turn,  were  employed  to  estimate  profile  drag.  Drag  is  the 
performance  quantity  that  is  most  sensitive  to  the  details  of  the  pressure 
distribution.  It  must  be  admitted  that  the  present  attack  is  speculative  in 
the  sense  that  the  data  available  are  limited  and  that  one  of  the  items  of 
great  interest,  the  drag,  was  not  measured  dxiectly  in  the  rotor  blade  tests. 

There  are  limitations  on  the  information  that  can  be  derived 
from  the  measured  rotor  pressure  differentials.  For  example,  the 
Reynolds  numbers  and  Mach  numbers  are  not  independent,  these  being 
related  to  the  radial  position  of  the  section  on  the  blade.  Certain  flow 
phenomena  are  practically  inseparable  with  respect  to  their  effects  on 
the  differential  pressure  distribution.  For  example,  the  induced  velocity 
distribution,  compressibility  effects,  angle  of  attack,  structural  distor¬ 
tion,  and  the  distortion  of  the  flow  about  the  blade  section  due  to  boundary- 
layer  growth  can  all  affect  one  or  more  of  the  pressure  distribution 
components  in  the  same  way  and  it  would,  indeed,  be  difficult  to  separate 
these  effects. 

The  discussions  that  follow  are  presented  in  the  context  of 
incompressible  flow.  However,  a  simple  Prandtl-Glauert  transformation 
is  introduced  where  necessary  to  correct  for  compressibility.  Such  a 
correction  (for  the  static  case  in  which  the  geometry  of  the  surface  is 
fixed)  merely  scales  all  the  terms  of  the  pressure  distribution  expansion 

by 

The  accuracy  of  the  experimental  data,  and  the  limitation 
imposed  by  the  number  of  points  at  which  the  differential  pressures  were 
measured,  influence  the  accuracy  that  can  be  obtained  in  the  decomposi¬ 
tion  of  the  pressure  differentials.  Care  has  been  exercised  in  an  attempt 
to  avoid  the  assignment  of  significance  to  fictitious  pressure  distribution 
components.  An  error  analysis  was  made  to  determine  at  what  point  the 
experimental  errors  can  introduce  a  distorted  view  of  the  components  of 
the  pressure  distribution. 

The  above  reservations  do  not  invalidate  the  present  approach. 
They  have  been  introduced  since  many  of  these  factors  reflect  the  diffi¬ 
culties  inherent  in  obtaining  the  measurement  per  se  rather  than  a  lack 
of  understanding  of  the  underlying  physical  phenomena.  Some  aspects, 
however,  are  still  beyond  the  state  of  the  art;  e.  g. ,  treatment  of  dynamic 
stall  and  compressibility  effects  in  regions  of  strong  spanwise  pressure 
gradients. 


2 


STATIC  ANALYSIS  OF  PRESSURE  DISTRIBUTIONS 


i 


THEORETICAL  TWO-DIMENSIONAL  PRESSURE  DISTRIBUTIONS 


In  linearized  two-dimensional  airfoil  theory,  the  pressure 
distributions  due  to  thickness,  angle  of  attack,  and  camber  can  be  calcu¬ 
lated  separately.  The  lift  and  the  moment  arise  only  from  angle  of  attack 
and  camber  because  the  pressures  on  the  upper  and  lower  surfaces  due  to 
the  thickness  cancel.  The  relationships  between  the  circulation,  lift, 
moment,  angle  of  attack,  and  camber  effects  will  be  expressed  in  terms 
of  so-called  "Glauert  Series"  (Reference  3). 

Let  the  airfoil  reference  chord  be  located  on  an  axis  designated 
the  %-a.xis,  with  the  leading  edge  at  -b  and  the  trailing  edge  at+b,  where 
b  is  the  half  chord.  Let  X  =  -b  cos  &  (Figure  2);  then  the  local  bound 
vorticity  strength,  ,  can  be  expressed  in  the  Glauert  Series  as 
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oo 

(t)  cot  ~  +21 


An(t)  sin 
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where,  in  general,  the  coefficients  are  functions  of  time,  and  the  steady 
two-dimensional  airfoil  problem  is  a  special  case.  The  total  circulation 
on  the  airfoil  is  given  by 


Thus, 


r 


d  X 


(2) 


r  =  2ir  o  (flo  +  l-R') 


(3) 


The  chordwise  distribution  cf  the  pressure  difference  is  given  by  the 
following  form  of  the  linearized  Bernoulli  equation  for  unsteady  flow; 
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J  A  jo(x,t) 
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(4) 
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where  ^  =  air  density,  V,  =  free-stream  velocity  at  infinity.  The  pres¬ 
sure  differential  can  be  expressed  in  terms  of  the  Glauert  coefficients  by 
substituting  (i)  into  (4).  This  result  is 


t)  =  2  p  Vt  I  R„  cot  A*.  sm  *1  & 

L  ^  v’-*' 

+  2  ?  b  *  7  A,)©  +  (R.+  j  At)  sm  0 

+  j  £  n  +  54,1  n9 

nst 


(5) 


The  airfoil  unsteady  lift  and  pitching  moment  about  the  midchord  per  unit 
span  can  also  be  expressed  in  terms  of  Glauert  coefficients  by  the  substi¬ 
tution  of  (5)  into  the  following  expressions: 


Thus, 


/ 

£_ 


«)  dX 


M 


b 

X  X)  dX 


L  ■  27Ti,ev,  [*.♦* A, +  ^r£(3/l.*  #,♦£().)' 
M  =  n?b‘V,  [fl.fi  A,-  A-  A  (A.fffl,- i-A,) 


(6) 

(7) 


It  is  observed  that  the  lift  depends  only  on  R0  and  R,  for  the 
steady  case  and  on  A,  ,  A,  ,  and  Rz  for  the  unsteady  case.  The  moment 
depends  on  Ra  and  Rz  for  the  steady  case  and  R0  ,  A,  ,  A,,  and  A  j 
for  the  unsteady  case.  The  physical  meaning  of  and  A,  in  thin  air¬ 
foil  theory  is  that  the  lift  contribution  from  the  A,  term  is  equivalent  to 
the  lift  of  a  flat  plate  at  a  given  angle  of  attack,  and  the  A,  term  is  equiva¬ 
lent  to  the  lift  due  to  camber  where  the  camber  line  chord  is  at  zero  angle 
of  attack.  ^  corresponds  to  the  pitching  m<  ment  at  zero  angle  of  attack. 


It  has  been  shown  in  Reference  3  that  the  steady  Glauert 
coefficients  R ^  are  related  to  the  camber  line  function  and  angle 


4 


of  attack  cL  ,  i.  e. , 


The  effect  of  the  presence  of  any  deflected  trailing  edge  tabs  (Figure  2) 
has  also  been  shown  to  be  (Reference  3) 


<  '  C,  *  2  %(  7  +  Stn  q  ) 


(10) 


where 


=  angle  of  deflected  trailing  edge  tabs 
(positive  downward) 

=  the  transformed  angular  position  through 
which  the  tab  is  rotated. 

Equation  (10)  may  be  rearranged  as 


(11) 


To  take  into  account  the  compressibility  of  the  air,  a  simple  Prandtl- 
Glauert  transformation  can  be  introduced.  Such  a  correction,  for  the 
steady  case  in  which  the  geometry  of  the  surface  is  fixed,  merely  scales 
all  the  terms  of  the  pressure  distribution  by  7/V?  -  M*  ’  where  is  the 
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Mach  number  in  the  free-stream.  This  transformation  will  not  hold  for 
the  unsteady  case. 

The  linearized  theory  predicts  an  infinite  pressure  differential 
at  the  leading  edge  of  a  flat  plate  of  zero  thickness.  It  can  be  seen  in 
Equation  (5)  that  as  &  approaches  zero  (the  leading  edge),  the  first 
term  in  the  Glauert  series  approaches  infinity.  This  reflects  the  inability 
of  the  linearized  theory  to  describe  the  large  flow  angle  change  (almost 
360°)  at  the  leading  edge  of  a  sharp  flat  plate  (leading  edge  radius  =  0). 

The  NACA  0012  and  0015  airfoils  have  finite  leading  edge  radii  and  finite 
pressure  differentials. 

Based  on  the  above  discussions,  the  steady  Glauert  coefficients 
for  the  NACA  0012  (Reference  4)  and  0015  (Reference  5)  airfoils  may  be 
computed  by  using  Equations  (8),  (9),  and  (11).  It  is  noticed  that  the 
NACA  0012  and  0015  airfoils  are  symmetric  and,  therefore,  the  camber 
line  function  =  0.  Table  I  shows  the  leading  coefficient  in  the 
Glauert  series  for  the  NACA  0012  airfoil  at  various  angles  of  attack. 

Table  II  shows  both  and  ft,  for  the  same  airfoil  as  Table  I  but  with 
3  degrees  of  flap  deflection  upward  at  trailing  edge.  Then's  for  the  NACA 
0015  airfoil  at  different  angles  of  attack  are  tabulated  in  Table  III.  The 
effects  of  compressibility  and  finite  leading  edge  radius  are  included  in 
all  the  tables.  The  analyses  in  References  6  and  7  were  used  to  calcu¬ 
late  the  effects  of  thickness,  finite  leading  edge  radius,  and  finite  trailing 
edge  angle.  Zero  angle  of  attack  data  are  snown  on  Figures  3  and  4  and 
are  discussed  in  a  later  section. 


MEASURED  TWO-DIMENSIONAL  PRESSURE  DIFFERENTIALS  (WIND- 
TUNNEL  DATA  FOR  NACA  0012  AND  NACA  0015  AIRFOILS) 

The  two-dimensional  pressure  differential  measurements  for 
NACA  0012  and  0015  airfoils  in  a  wind  tunnel  are  reported  in  References 
4  and  5,  respectively.  The  ranges  of  angle  of  attack  and  Mach  numbers 
are  indicated  in  Tables  I,  II,  and  III.  Locations  (in  percent  chord)  of  the 
pressure  taps  aft  of  the  leading  edge  are  tabulated  below. 


NACA  0012  AIRFOIL  (REFERENCE  4) 


TAP  NO.  1  2  3  4  5  j 

6  j 

7  j  8  1  9 

10 

II 

12 

13 

14 

—  -  -  .  v  +-  t  1 

1  j 

POSITION  1  ,  „„„ 

x/c  0. 008 j 0.017] 0.040 j  0.065 jO. 090 

jO.  130 

0. 108  0.233  0.335 

1  1  1 

0.500 

0.625 

0.769 

0.915 

0.960 

NACA  0015  AIRFOIL  (REFERENCE  5) 
TAP  NO.  '  1  2  3  4  I 


5  j  6 


POS'!!°N'0. 025^  0. 05.0. 10*0.  i5!0.20!o.25jo.  30|0.35l0.  '10 

x'  4  |  !  i  1  I  i  I  1  I 


10 

j  11 

12 

13 

L"> . 

15  1 

l 

Ji! 

0.45 

0.50 

1 

0.55 

0.60 

0.70 

0.80 

0 . 90  j 
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Extraction  of  the  Glauert  coefficients  from  the  measurements;  by  means 
of  Equation  (5)  is  subject  to  limitations  such  as  the  following: 


1.  With  a  finite  number  of  pressure  taps  along  the 
chord  of  the  airfoil  section,  it  is  impossible  to 
calculate  all  the  Glauert  coefficients  in  Equation 
(1).  The  series  must  be  truncated  at  some  finite 
number,  W  ,  less  than  or  equal  to  the  number  of 
taps. 

2.  The  accuracy  of  the  experimental  results  and 
the  numerical  errors  in  the  procedure  for  com¬ 
puting  the  Glauert  coefficients  should  be  considered. 


It  may  be  reasonable  to  determine  the  accuracy  of  the  reading 
at  each  pressure  tap  or  the  possible  percentage  error  of  that  tap.  Then, 
the  number,  N  ,  of  terms  to  be  retained  in  the  series  can  be  determined 
such  that  the  contribution  due  to  thefN  +  j)tlu  term  in  the  Glauert  series 
at  a  given  pressure  tap  is  within  the  probable  error.  In  other  words,  the 
upper  and  lower  bound  can  be  established  for  each  pressure  tap  depending 
on  the  possible  percentage  error.  From  the  experimental  results,  a 
pressure  differential  curve  can  be  drawn  from  the  readings  of  the  pressure 
taps  along  with  the  upper  and  lower  bound?  of  this  pressure  differential 
curve  (see  Figure  5).  The  number  N  will  be  determined  in  such  a  fashion 
that  the  first  W  terms  of  the  Glauert  series  lie  between  the  curves  of  the 
uppe;  and  lower  bounds. 

Unfortunately,  the  problem  of  accuracy  for  each  pressure  tap 
has  not  been  discussed  in  References  4  and  5.  However,  the  pressure 
differentials  at  each  chordwise  location  measured  at  zero  angle  of  attack 
will  provide  useful  information  on  the  accuracy;  for  example,  see  the 
experimental  data  given  by  the  first  table  on  page  62  of  Reference  5. 
Theoretically,  if  the  pressure  tap  measurements  on  the  upper  and  lower 
surfaces  are  all  perfect,  the  measured  pressure  differentials  should  be 
equal  to  zero  at  zero  angle  of  attack  since  the  NACA  0015  airfoil  is 
symmetric.  The  values  of  the  pressure  differences  provide  an  estimate 
of  the  experimental  accuracy.  The  concept  of  a  so-called  "root  mean 
square  value"  may  be  borrowed  from  numerical  analysis  to  define  the 
accuracy  instead  of  the  possible  percentage  error.  The  concept  of 
possible  percentage  error  will  be  discussed  in  detail  in  the  next  section. 

The  root  mean  square  value  is  defined  as 


£  = 


rms 
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where  A^s-  is  the  pressure  differential  measured  at  zero  angle  of 
attack^and  K  is  the  total  number  of  differential  pressure  taps. 

The  root  mean  square  values  for  References  4  and  5  at  different 
Mach  numbers  are  tabulated  respectively  as  follows: 


NACA  00' 2  AIRFOIL  (REFERENCE  4) 

MACH  NO. |  0.30  jO.IO  jO.SO  'o.6Q  0.65  i 0.70  j  0.75  jp.80  ; 

e  |  0.054  0.055]  0.0681  0.039!0.04li0.040l  0.036  0.053 
rros  .  I  i  >  i  t 


NACA  0015  AIRFOIL  (REFERENCE  5) 
MACH  NO. | 
f  rms 


1 0. 6S0{  0 . 675 

I 

0.701 

0.751 

0.777 

0.805 

0.831 

0.040j0.042 

0.018 

0.058 

0.063 

0.096 

0.090 

The  at  Mach  numbers  other  than  those  listed  in  the  above  table  may 

be  estimated  by  interpolation. 

Once  6rrn%  is  given,  it  is  required  that  the  error  due  to  numeri¬ 
cal  computation  should  be  smaller  than  the  given  .  In  other  words. 


rm  i 


^  e 

rms 


where  is  the  pressure  differential  by  the  Glauert  series  (see 
Equation  5). 

The  method  selected  for  computing  the  Glauert  coefficients  was 
the  "least-square  polynomial  approximation"  This  method  has  the  advan¬ 
tage  of  using  all  the  given  information  from  the  pressure  taps  and  of 
providing  the  required  results  with  a  minimum  effort  in  computation.  The 
details  of  this  method  are  given  in  Appendix  I.  The  following  example  is 
given  to  demonstrate  the  method. 
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Airfoil  type  NACA  0015 

Mach  number  0.  701 


Angie  of  attack  -2  degrees 


Let  n  be  the  number  of  terms  to  be  calculated  in  the  Glauert 
series  by  the  method  of  "least-square  polynomial  approximation".  One 
obtains 


n 

5 

6 

7 

8 

9 

erms 

0.070 

0.057 

0.037 

0.036 

0.031 

to  compare  with  6rms  =  0.0476  {based  on  the  pressure  differentials  at  zero 
angle  of  attack  and  Mach  number  0.701).  The  implication  of  this  example 
is  that  if  all  the  information  of  the  pressure  differentials  is  used  to 
compute  the  first  five  coefficients  in  the  Glauert  series,  the  root-mean- 
square  error  between  and  A^ris  larger  than  the  estimated  root-mean 

square  error  due  to  the  experimental  errors.  When  the  first  seven  terms 
are  used,  however,  the  £'  *s  sma^er  than  6rms  as  expected.  Then, 
the  computation  may  be  stopped  at  k  =  7  for  the  present  case,  since  no 
useful  information  will  be  provided  by  the  extra  terms.  It  should  be 
pointed  out  that  the  value  of  Yl  depends  on  £r„s  and  will  be  different 
for  various  £r«,*  • 

The  results  of  the  numerical  computation  for  NACA  0012  and 
0015  airfoils  are  presented  in  Tables  IV  through  VIII.  Only  R0  ,  R(  ,  and 
are  presented,  since  lift  and  moment  depend  only  on  those  coefficients 
in  steady-state  condition. 

MEASURED  ROTOR  PRESSURE  DISTRIBUTIONS  (FLIGHT  DATA  FOR 
H-34  WITH  NACA  0012  AIRFOIL  AND  UH-1A  WITH  NACA  0015  AIRFOIL) 


The  coefficients  of  the  series  representations  of  the  measured 
differential  pressure  distributions  for  the  H-34  rotor  blade  with  NACA 
0012  airfoil  (References  8  and  9)  and  the  UH-1A  rotor  blade  with  NACA 
0015  airfoil  (Reference  10)  have  been  computed  for  the  following  cases: 
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ADVANCE  RATIO 

RADI ALSTATIOK 
r 


H-34 


0.18,  0.29 


0.40,  C.8S,  0.90 


UH-1A 


0.0,  0.08,  0.26 


0.40,  0.85,  0.90 


(H-34)  0.042  0.158  0.300  0.600  0.910 
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where  the  summation  is  truncated  according  to  the  number  N  of 
pressure  differential  measurement  stations, and  the  subscript  JL  and  -k. 
indicate  the  i“'  and  ^“position  of  r  and  Iff  ,  respectively. 

The  evaluation  of  the  accuracy  in  reading  the  oscillograph 
records  or  the  pressure  transducers  has  been  discussed  in  References  8, 
9, and  10, respectively.  In  References  8  and  9,  the  estimated  accuracy  of 
each  individual  data  point  of  the  poorest  records  is  +  5  percent  with  99.  7 
percent  confidence,  i.  e. ,  99.  7  percent  of  points  have  an  error  less  than 
+  5  percent.  It  is  stated  in  Reference  10  that  overall  errors  in  the  tabu¬ 
lated  data  of  +  5  to  +  7  percent  are  to  be  expected,  but  a  set  of  figures  for 
pressure  transducer"  calibration  curves  is  also  given  in  Reference  10.  The 
pressure  differences  between  measured  and  computed  data  may  then  be 
plotted  to  indicate  the  accuracy  of  the  computed  data. 

When  the  method  of  "least-square  polynomial  approximation"  is 
applied  to  compute  the  Glauert  coefficients  in  Equation  (12),  the  possible 
percentage  errors  of  +  5  percent  are  required  for  Preferences  8  and  9, 
while  +  5  to  +  7  percent  are  required  for  Reference  10.  An  example  for 
each  case  is~given  for  demonstration. 

t- 

The  first  case  to  be  considered  is  the  pressure  differentials 
measured  at  /x  =0.  18,  r  =  0.  90,and  Iff  =  0  degrees  (Reference  8). 


x/c 

4P(lb/in?) 

A?'n  «  a 
n  *  3 

ePP,  (« 

APn  <=  4 

€PP,  00 

AP' 
n  =5 

€ppe  (50 

0.017 

7.715 

7.786 

0.93 

7.744 

0.13 

7.720 

0.68 

0.090 

3.981 

3.739 

6.09 

3.843 

3.47 

3.925 

1.15 

0.168 

2.858 

2.869 

3.77 

2.934 

2.65 

2.950 

3.22  | 

0.233 

2.455 

2.441 

0.58 

2.  4^3 

0.49 

2.397 

2.35  | 

0.335 

1.741 

1.920 

10.30 

1.826 

4.89 

1.744 

1.49  ! 

0.625 

0.635 

0.6425 

1.17 

0.542 

0.473 

0.639 

0.55 

0.015 

0.042 

-0.515 

222.61 

0.103 

155.24 

0.0404 

3.87 

where  A^5n>j=  computed  pressure  differentials  by  using  first  Yl  (=  3)  terms 
of  Glauert  series,  and 


possible  percentage  error 
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From  this  table. 

inside  the  upper  and  lower  b  Glauert  series  give  a  repre¬ 

sentational  the  ^nea  sur^d  datiScons  is  t  ent  with  enperitnental  accuracy. 

?'rerrshow^the  difference  between  c-putedjressure^d^re^ 
KSurS  ^MTi-5.  «»  curves  given  by  ±  7  percent. 

The  rest  of  the  tneasured  pressures  ar^treated^th^same 

fashion  to  yield  the  first  n  terms^  f  h  ^  NACA  0012  and  0015  airfoils 
of  C0l(Ca+  ?C,)  +  £  ,)  lenprtivelv  The  purpose  of  using 
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THEORETICAL  UNSTEADY  ROTOR  PRESSURE  DISTRIBUTIONS 


The  theoretical  unsteady  Glauert  coefficients  for  the  flight 
conditions  to  be  considered  were  calculated  in  Reference  1.  The  flight 
conditions  are:  H-34  rotor  blade  (NACA  0012  airfoil)  at  p-  =  0.  18  and 
p.  =  0.29,  and  UH-lA  rotor  blade  (NACA  0015  airfoil)  at  p.  =  0.08  and 
p.  =  0.  2b.  Although  details  of  the  calculations  are  given  in  the  cited 
reference,  the  basic  assumptions  are  repeated  here  to  assist  in  the 
interpretation  of  the  computed  unsteady  Glauert  coefficients.  They  are 

1.  The  wake  configuration  can  be  adequately  prescribed. 

2.  The  spanwise  blade  slopes  (e.  g.,  coning,  bending),  and  the 
section  angles  of  attack  below  stall  are  small. 

3.  The  in-plane  components  of  the  induced  velocities  at  the 
tip-path  plane  are  small  and  can  be  neglected. 

4.  Belov/  stall,  the  lift-curve  slope  is  constant. 

5.  The  blade  section  circulation. is  limited  to  a  maximum 
value  for  angles  of  attack  at  and  above  stall. 

6.  For  angles  of  attack  above  stall,  the  blade  section's 
normal  force  is  taken  as  the  sum  of  the  stall-limited 
circulatory  force  and  a  cross-flow  drag  force. 

7.  The  Mach  number  and  Reynolds  number  effects  are 
assumed  to  influence  only  the  lift-curve  slope. 

8.  The  interference  effects  of  the  rotor  hub,  fuselage, 
etc. ,  are  negligible. 


The  pressure  distributions  and  Glauert  coefficients  are 
related  through  Equation  (5).  The  computed  R^s  are  presented  in 
Figures  7  through  18. 

MEASURED  ROTOR  PRESSURE  DISTRIBUTIONS 

The  Glauert  coefficients,  Cj  ,  for  the  measured  rotor  pressure 
differentials  have  been  calculated  at  a  given  set  of  radial  and  azimuthal 
stations  from  Equation  (12).  The  C^'s  cannot  be  compared  directly  with 
the  obtained  from  Reference  1.  However,  the  relationships  between 

the  Cj-s  and  s  can  be  established.  First,  it  is  noted  that  the  rate  of 
change  of  azimuthal  station  is  the  shaft  speed;  i.  e.  , 


13 


(13) 


iJt  =  oJ 
d  t 

where  tJ  =  rotational  speed  of  the  rotor  blade. 


the  identity 


If  Equation  (13)  is  substituted  into  Equation  (5)  with  the  aid  of 


&  -  2  (  sin  9  -  j  sin  2  &  *  y  5in  3  9  - - ) 


-  7T  <  9  <  rr 


Equation  (5)  reduces  to 


A  p(&yyr)  -  2  p  V;  (Rocot-Y  +  LRn5inYL&) 
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n  +  «  (14) 
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If  the  coefficients  of  like  trigonometric  functions  are  collected, 
the  following  set  of  equations  that  relate  the  CA's  of  Equation  (12)  to  the 
ftx's  of  Equation  (14)  is  obtained  [by  assuming  Af>(8,t)  »  A f>(6,  y^)  for 
periodic  airloads]. 


C,iyr )  -  R,(p  =■ 


Cjyr)-ftXp  « 


fljyr) 
b 


?  +  I*  Sin  yr  dy 


3  flip) +Q,a/r)  +  -jr 


(15) 
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b  = 
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Once  the  £  '5  are  given,  then  the  's  can  be  found  by  solving 
the  set  of  simultaneous  differential  equations  (Equation  (15)).  The  solu¬ 
tion  to  Equation  (15)  is  straightforward  if  both  Cz 's  and  '*  are  expanded 
in  Fourier  series  with  respect  to  azimuthal  angle  p"  ,  i.  e.  , 


/*t  *  /£  +  2  n  ^ 

H.  r  t 


+  Co?  n  Tfr) 


+  cos  n 


(16) 


The  number  of  terms  H.  has  to  be  truncated  at  some  finite  limit  for 
practical  computational  considerations.  In  the  present  analysis,  =  5  is 
used  since  the  lifts  and  moments  are  dependent  on  the  first  four  Glauert 
coefficients  (  A0  ,  A,  ,  fit,a.nd  )^  The  process  of  substituting  Equation 
(16)  into  (15)  and  equating  the  coefficients  of  like  trigonometric  functions 
again  yields  a  set  of  algebraic  equations  (see  Appendix  II)  which  can  be 
solved  readily  on  a  digital  computer."'  The  results  for  the  first  three  /Vs 
are  plotted  in  Figures  7  through  18. 


I  4 


*The  Fourier  expansion  for  Cj%  in  ^  is  a  standard  subroutine  in 
programming. 
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COMPARISON  AND  DISCUSSION  OF 
PRESSURE  DISTRIBUTION  COMPONENTS 


It  must  be  recalled  for  purposes  of  this  discussion  that  a 
particular  method  of  curve  fitting  was  used;  namely,  least-square  fit. 
Furthermore,  it  must  be  noted  that  the  function  representing  the  pressure 
distributions  were  not  all  orthogonal- -the  angle-of-attack  effect  being 
contained  in  a  term  coif  t  which  is  not  orthogonal  with  respect  to  the 
rest  of  the  Glauert  series  terms.  The  important  aspect  of  this  observation 
is  that  a  change  of  the  number  of  points  to  be  fit  or  of  the  number  of  terms 
to  be  retained  can  produce  a  new  set  of  coefficients? 

The  Glauert  coefficient  comparisons  that  follow  will  be  grouped 
in  accordance  with  this  order; 

•  Two-dimensional  theoretical  compared  to  two- 
dimensional  measured  coefficients;  effects  of  viscosity 

•  H-34  pressure  distributions  as  functions  of 
spanwise  station 

a  Comparison  of  H-34  data  at  an  advance  ratio 
of  /x  =  0. 18  with  data  at  n-  -  0.  29 

•  UH-1A  pressure  distributions  as  functions  of 
spanwise  station 

•  Comparison  of  UH-1A  data  at  advance  ratios  of 
/x  =0.08  with  data  at  yo  =0.  26 

•  Comparison  of  H-34  data  with  UH-1A  data 

•  Comparison  of  two-dimensional  data  with  rotor  blade  data 

TWO-DIMENSIONAL  RESULTS;  EFFECTS  OF  VISCOSITY 

Theoretical  Glauert  coefficients  based  on  linearized  analysis 
for  an  inviscid  fluid  (Reference  3)  are  presented  in  Tables  I,  II, and  III. 
These  are  the  coefficients  representing  the  pressure  differentials  for  the 
steady  flow  case  for  the  following  configurations: 

Table  I  -  NACA  0012,  0&<  or  «  16°  ,  0.  3 * /%&  0.  8 

Table  II-  NACA  0012,  3°  Negative  Trailing  Edge  Flap, 

0%^  *  16:  0.  3  #  m  <  0.  8 

1  '00 

Table  HI-  NACA  0015,  2°  -  <*  <  12°  ,  0.  3  «  ff,*  0.84 
*Each  coefficient  of  an  orthogonal  series  can  be  determined  independently. 
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Calculations  contained  in  Table  II  were  carried  out  for  the  flap  configura¬ 
tion  used  on  the  outboard  sections  of  the  H-34  blades  (Reference  9)  and  in 
the  two-dimensional  wind-tunnel  tests  reported  in  Reference  4. 

Within  the  limits  of  the  linearized  theory  used,  the  pressure 
differentials  on  the  NACA  0012  and  0015  sections  are  the  same  (for  the 
same  angle  of  attack,  ,  and  Mach  number,  )  because  they  both  have 
zero  camber.  The  corresponding  theoretical  pressure  differentials  have 
only  Ra  amplitudes  (see  Tables  I  and  III).  A  deflected  flap  introduces  an 
effective  camber  and,  consequently,  an  f 1 ,  pressure  differential  compo¬ 
nent  and  a  shift  in  the  angle  of  attack  for  zero  lift.  This  accounts  for  the 
differences  between  Tables  I  and  II. 

Results  of  least-squares  fits  to  the  pressure  differential  data 
given  in  References  4  and  5  are  shown  in  Tables  IV  through  VIII.  Compari¬ 
sons  of  Table  I  with  Table  IV,  Table  II  with  Table  VI,  and  Table  III  with 
Table  VIII  immediately  establish  the  lack  of  agreement  between  the  theo¬ 
retical  and  experimental  values.  One  of  the  more  striking  differences  is 
the  presence  of  large  Rz  components  in  the  experimental  data. 

The  Ri  coefficient  data  of  Tables  IV  through  VIII  were  normalized 
with  respect  to  the  sum  (  R„  *■  jr  R,  )  since  this  sum  is  proportional  to  the 
net  airfoil  circulation,  and  the  results  are  presented  in  Tables  IV  through 
VIII  as  normalized  quantities,  Ri  .  The  lack  oi  consistency  noted  in  the 
R L  values  leads  to  the  tentative  conclusion  that  important  viscous  and 
compressibility  effects  may  be  present  and,  possibly,  interactions  between 
these.  No  precise  boundary  separates  the  two  effects. 

It  is  well  known  that  the  viscous  effect  for  a  flow  over  a  smooth 
body  will  be  confined  to  a  very  thin  layer,  the  so-called  boundary  layer. 

The  flow  outside  the  boundary  layer  is  a  potential  flow  over  a  slightly 
different  body  which  consists  of  the  original  airfoil  and  an  additional 
thickness  distribution  due  to  the  viscous  effect.  The  NACA  0012  and  NACA 
0015  airfoils  are  symmetric;  i.  e. ,  they  have  no  geometric  camber. 
However,  viscous  effects  will  modify  the  symmetry  of  those  airfoils. 

Thus,  a  camber  line  due  to  the  viscous  effect  is  induced  and  contributes 
to  the  value  of  the  R,  and  higher  Glauert  terms  observed  in  the  two- 
dimensional  wind-tunnel  data.  As  speed  increases,  the  boundary  layer 
may  change  from  laminar  to  turbulent.  If  the  flow  becomes  turbulent, 
boundary  layer  displacement  will  be  much  thicker  than  the  laminar  one. 

If  the  flow  separates,  the  boundary  layer  displacement  will  also  be 
thicker.  Both  cases  will  increase  the  magnitude  of  A,  . 

An  increase  in  free-stream  velocity  increases  the  local  Mach 
number  for  the  profile  until  a  zone  of  local  supersonic  flow  has  appeared 
and  shock  waves  occur.  As  a  result  of  the  very  large  adverse  pressure 
gradients  imposed  by  the  shock  wave,  the  boundary  layer  thickens  rapidly, 
generally  becomes  turbulent,  and  often  separates  from  the  profile,  thus 
producing  a  large  increase  in  the  thickness  of  the  wake.  The  change  in 
effective  body  shape  usually  acts  to  decrease  the  lift  and  increase  the  drag. 
This  coupling  of  the  interaction  between  the  shock  wave  and  boundary  layer 
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presents  a  formidable  problem  (Reference  11). 

Boundary  layer  separation  may  develop  at  the  leading  edge  or 
progress  forward  from  the  trailing  edge.  Leading  edge  separation  will  be 
emphasized  here.  Types  of  leading  edge  separation  occurring  on  airfoils 
may  be  classified  roughly  as  (1)  separation  with  laminar  reattachment,  (2) 
separation  with  turbulent  reattachment,  and  (3)  separation  without  reattach¬ 
ment.  As  far  as  NACA  0012  and  NACA  0015  airfoils  are  concerned,  the 
most  likely  type  in  the  operating  lift  coefficient  range  is  type  (2).  The 
flow  over  a  sharply  curved  leading  edge  encounters  a  strong  adverse  pres¬ 
sure  gradient  at  high  angle  of  attack.  The  flow  cannot  travel  around  the 
bend  without  laminar  separation  taking  place.  After  separation,  transition 
to  turbulence  occurs  in  the  boundary  layer  somewhere  downstream  of 
separation.  Transport  mechanisms  associated  with  the  turbulent  layer 
may  cause  reattachment.  In  between  the  separation  and  attachment  points, 
there  is  an  enclosed  bubble  of  air.  The  size  of  the  bubble  has  been  found 
(Reference  12)  to  vary  from  a  very  small  fraction  of  the  chord  to  some¬ 
thing  comparable  with  the  chord  length.  Airfoil  thickness  was  found  to  be 
a  parameter.  The  tests  in  Reference  12  show  that  a  thin  wing  at  small 
angle  of  attack  develops  a  very  small  bubble  near  the  leading  edge.  The 
size  of  the  bubble  increases  with  angle  of  attack  until  it  encompasses  the 
whole  upper  surface.  For  the  moderately  thicker  wing,  the  bubble  con¬ 
tracted  slightly  up  to  a  certain  angle  of  attack;  beyond  this,  the  bubble 
suddenly  burst  and  caused  an  abrupt  stall. 

The  presence  of  the  air  bubble  and  the  boundary  layer  displace¬ 
ment  will  change  the  effective  shape  of  the  airfoil  profile.  Size  of  the  air 
bubble  and  the  thickness  of  boundary  layer  displacement  depend  on  the 
speed  of  the  freestream.  With  these  two  points  in  mind,  it  is  not  sur¬ 
prising  that  in  Tables  IV  and  VIII  the  ft,  terms  are  dominant.  Theoretically, 
f\a  increases  as  both  angle  of  attack  and  Mach  number  increase.  However, 
the  results  of  using  a  least- square  fit  to  the  experimental  data  show  that 
the  term  of  the  NACA  0012  airfoil  increases  with  angle  of  attack  but 
decreases  with  Mach  number  (see  Table  IV).  The  term  of  NACA  0015 
airfoil  is  quite  random  with  respect  to  angle  of  attack  and  Mach  number. 

No  easily  discernable  patterns  are  apparent  in  Table  VIII. 

Lack  of  trends  in  the  two-dimensional  pressure  distributions 
will  make 'the  rotor  blade  data  obtained  in  flight  difficult  to  interpret.  It 
is  worthwhile,  therefore,  to  delve  further  into  sources  of  the  trouble.  It 
appears  that  the  source  of  the  scatter  might  be  the  method  of  testing,  i.  e. , 
changes  of  speed  introduce  both  Reynolds  number  and  Mach  number  changes. 
Further,  no  wall  corrections  were  made  for  the  0015  pressure  data. 

The  effects  of  change  in  Reynolds  number  and/or  angle  of 
attack  on  the  separation  characteristics  of  a  two-dimensional  airfoil 
can  be  estimated  with  the  aid  of  a  "Ville  Plot"  (Reference  13)  such  as 
that  shown  in  Figure  19. 
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Figure  19  shows  the  influence  on  leading  edge  separation  of 
the  nondimensional  leading  edge  radius,  j8  ,  angle  of  attack,  «c  ,  and 
Reynolds  number  based  on  leading  edge  radius,  R,  .  The  effects  of 
compressibility  are  not  displayed.  The  importance  of  this  plot  with 
regard  to  the  Glauert  coefficients  is  that  the  effective  shape  and,  hence, 
the  coefficients  are  changed  as  the  operating  point  (speed,  angle  of 
attack)  is  changed. 

In  essence,  the  study  of  the  two-dimensional  airfoil  data  by 
the  technique  of  developing  Glauert  coefficients,  by  least-squared  curve 
fitting,  has  revealed  several  unexpected  aspects  with  respect  to 
the  interpretation  of  the  data.  It  will  be  seen  that  a  similar  situation 
exists  for  the  rotor  data. 

H-34  DATA  AT  ADVANCE  RATIO  OF  p,  =0.18 


Figures  7,8,  and  9  show  the  first  three  nondimensionalized 
Glauert  coefficients  for  the  H-34  flying  at  an  advance  ratio  of  p.  =  0.  18 
at  radial  stations  r/R  =  0.40,  0.85,  and  0.  90,  respectively.  The  theo¬ 
retical  curves  identified  by  the  short  dashed  lines  were  obtained  from 
Reference  1  .  The  solid  line  marked  "experimental  data"  corresponds 
to  the  determination  of  the  Glauert  coefficients  from  pressure  distribu¬ 
tions  for  the  unsteady  case ,  that  is,  the  representation  (Equation  (5)). 
The  Ci  coefficients  (Equation  (12))  are  the  straightforward  fitting  of  the 
pressure  distributions  as  though  there  were  no  unsteady  effects. 

The  differences  between  the  three  curves  could  also  be  de¬ 
scribed  in  the  following  manner:  The  theoretical  data  are  based  on  the 
assumption  that  the  behavior  of  the  airfoil  sections  is  predominantly  that 
of  a  flat  plate  without  boundary-layer  growth,  whereas  the  experimental 
data  contained  the  viscous  effects  and,  of  course,  reflect  the  true  environ¬ 
ment  for  each  section.  The  difference  between  the  experimental  fit  for 
the  curve  and  the  curves  marked  Ct-  should  represent  the 

unsteady  effects. 

Figure  7  shows  the  theoretical  and  two  experimental  rep:  esenta- 
tions  of  the  nondimensional  pressure  distributions  for  the  H-34  at  advance 
ratio  of  p  =  0.18  and  at  radius  r/R  =  0.  40.  It  is  immediately  evident 
that  the  theoretical  predictions  would  indicate  very  small  values  for  all 
the  coefficients  except  .  This  is  not,  however,  borne  out  by  the 
experimental  data.  The  experimental  data  show  sizable  values  for 
on  the  retreating  side  of  tue  disc  and  sizable  values  of  over  virtually 
all  azimuthal  positions.  Furthermore,  there  appears  to  be  an  unsteady 
effect  on  Rc  and  ft,  on  the  retreating  side  of  the  disc;  that  is,  for  the 
region  180°  <  IF  <  360?  There  is  also  an  unsteady  effect  contributing 
to  around  the  complete  azimuth. 

Figures  8  and  9  show  the  azimuthal  distribution  of  the  coeffi¬ 
cients  at  nondimensional  radial  stations  r/R  =  0.  85  and  0.  90,  respectively, 
for  this  same  flight  condition.  Again,  the  experimental  data  do  not  agree 
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particularly  well  with  the  theoretical  data  for  R,  and  Rt ;  even  the 
flo  coefficients  for  the  experimental  data  depart  from  the  theoretical 
coefficient  values  on  the  advancing  side  of  the  disc. 

The  unsteady  effects  at  these  larger  radial  stations  (  r/R  = 

0.  85  and  0.  90)  appear  to  be  considerably  less  than  at  the  inboard  station 
r  /R=  0.4,  as  can  be  seen  by  comparing  the  differences  between  the 
experimental  Ri  coefficients  and  the  experimental  Cx-  coefficients  of 
Figures  8  and  9  and  those  of  Figure  7.  The  other  change  that 
appears  to  occur  with  radial  station  is  tnat  R,  and  RL  have  large  values 
on  the  advancing  side  at  the  larger  radial  stations  (Figures  8  and  9), 
whereas  these  coefficients  peak  on  the  retreating  side,  as  shown  on 
Figure  7,  at  the  inboard  station.  This  change  with  radial  station  may 
simply  be  a  manifestation  of  compressibility  and  tip  effects  dominating 
on  the  advancing  side  at  the  outer  radial  station  and  the  stall  effect  con¬ 
tributing  to  the  distortions  of  the  pressure  distribution  at  the  inboard 
radii. 

H-34  AT  ADVANCE  RATIO  fx  =  0.  29 

The  data  for  the  H-34  at  advance  ratio  of  /j-  =  0.  29  are  shown  in 
Figures  10,  11, and  12.  At  the  most  inboard  radial  station  for  this  flight 
condition  (Figure  10),  the  excursion  of  the  coefficients  on  the  retreating 
side  is  quite  large.  In  fact,  it  was  necessary  to  change  the  scale  of 
the  ordinate  on  Figure  10  compared  to  the  others  in  the  set  of  Figures  7 
through  12,  The  other  notable  aspect  about  Figure  10  is  the  occurrence 
of  quite  large  negative  values  for  and  C,  on  the  retreating  side  of 
the  disc.  The  coefficients  fl»  and  C.t  determined  from  the  experiment 
are  also  negative  for  part  of  the  retreating  side  of  the  disc.  Again  the 
outboard  sections  exhibit  approximately  the  same  characteristics  as 
those  observed  at  the  lower  advance  ratio;  namely,  unsteady  effects  _ 
are  relatively  small,  and  the  experimental  values  deduced  for  and 
are  considerably  larger  than  the  corresponding  analytic  coefficients. 


COMPARISON  OF  H-34  DATA  AT  AN  ADVANCE  RATIO  OF  R-  =0.  18 
WITH  DATA  AT  p.  =0.29 

Strong  similarities  exist  for  the  outboard  sections  of  the  H-34 
at  the  two  advance  ratios  investigated.  The  inboard  section,  on  the  other 
hand,  appears  to  have  notable  differences._  This  is  best  seen  on  the 
retreating  side  of  the  azimuth, where  the  A,  coefficients  are  of  different 
signs  in  the  two  cases  and  the  f\?  coefficients  also  differ  ia  sign  over 
a  portion  of  the  retreating  side.  Obviously,  the  proximity  of  the  blade 


All  three  measures  of  the  circulation- -that  is,  theoretical  ( fl0  +  jr  f%  ), 
experimental  (  Ro  +  ?  ft,  ),  and  flight  (  C„  *  %  C,  )--are  in  reasonable 
agreement  and,  hence,  the  azimuthal  variation  of  lift  is  similar.  This 
only  suggests  that  the  gross  flow  properties  (e,  g. ,  local  momentum)  are 
relatively  insensitive  to  viscous  and  unsteady  effects. 
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section  to  the  reverse  flow  region  could  introduce  the  change  in  observed 
behavior.  The  tangential  velocity  at  r/R  =  0.  4  differs  by  a  factor  of  two 
for  these  advance  ratios  at  =  270  degrees.  The  low  velocities  (about  70 
.ft/sec  for  fj.  =0.29  and  about  140  ft/sec  for  /*.  =  0.  18)  make  the  section 
characteristics  especially  sensitive  to  inflow  velocities  that  are  locally 
high.  Consequently,  the  local  angle  of  attack  and  the  local  variation 
of  the  velocity  over  the  chord  could  introduce  the  reversal  of  the  camber 
effect  noted  in  the  comparison  of  these  two  cases. 

UH-1A  PRESSURE  DISTRIBUTIONS  AS  FUNCTIONS  OF  SPANWISE 
STATION 

The  UH-1A  data  analyzed  are  shown  in  Figures  13,  14, and  15 
for  the  case  of  /j,  =  0,08  and  in  Figures  16,  17, and  18  for  =  0.26. 

_  _  Figure  13  displays  the  data  at  pu  =  0.08  and  r/R  =  0.4.  The 

ft0  ,  ft,  ,  ftt  ,  and  corresponding  C  are  relatively  constant  over  the 
azimuth  except  in  the  region  from  about  T/r  =150  degrees  to  about  y~  = 
230  degrees._  That  azimuthal  region  is  characterized  by_a  change  in  sign 
in  ft,  and  Ct  and  a  sizable  change  in  the  value  of  the  fta  coefficient. 

At  radial  stations  r/R  =_0.  85  and  0.  90  (Figures  14  and  15, 
respectively),  it  is  seen  that  the  fti  and  Ci  coefficients  representing  the 
experimental  data  lie  very  close  to  each  other  and,  hence,  indicate  that 
the  unsteady  effects  are  small.  The  other  item  of  interest  in  Figures  14 
and  15  is  that  the  experimental  ft0  coefficients  are  relatively  low  com¬ 
pared  to  the  theoretical  value  and  the  ft,  and_  ft2  coefficients  are  quite 
sizable,while  the  theoretical  coefficients  for  ft,  and  f\z  are  virtually 
zero.  Returning  again  to  Figure  13,  it  is  noted  that  the  agreement  between 
the  theoretical  and  experimental  data  is  reasonable  for  fto  but  has  appre¬ 
ciable  errors  in  the  ft,  and  ftt  values.  Again,  the  data  at  r/R  =0.4 
indicate  larger  unsteady  effects  than  do  the  data  for  the  outboard  stations. 

The  UH-1A  data  at  ju*  =  0.  26  are  shown  in  Figures  16,  17,  and 
18  at  radial  stations  r/R  =  0.  40,  0.  85, and  0.  90,  respectively.  Figure  16 
shows  that  the  Aa  component  dominates  on  the  advancing  side_of  the 
blade;  but  on  the  retreating  side  (180°  <  y  <  360^),  both  the  ft, 
and  ftz  coefficients  become  sizable.  Here,  again,  the  difference  between 
the  fy  and  Ci  coefficients  indicates  a  substantial  effect  traceable  to 
unsteady  aerodynamics.  Figures  17  and  18  indicate  sizable  excursions  in 
all  the  measured  coefficients  at  the  outboard  stations.  Even  the  fta  terms 
vary  over  a  range  from  about  0.  4  to  a  value  somewhat  greater  than  1. 0. 
Therefore,  the  theoretical  curves  which  are  dominated  by  fta  are  in 
disagreement  with  the  experimental  pressure  distribution  components. 

The  fl0  and  Pi,  coefficients  generally  behave  the  same  at  r/R  =,0.85 
and  0.  90,  and  similarly  the  coefficients  are  of  the  same  type." 


Note  the  change  in  the  scale  on  Figure  17. 
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COMPARISON  OF  UH-1A  DATA  AT  ADVANCE  RATIO  jn.  =  0.  08  WITH 
DATA  AT  ><.  =  0.  26  - - - 


Data  obtained  at  an  advance  ratio  of  M  =o-  26  seem  to  have 
larger  excursions  than  the  data  obtained  at  /j~  =  0.08  for  the  UH-1A  rotor. 
The  shapes  of  the  respective  curves  at  r/R  =  0..  85  and  0.90  are  roughly 
similar, but  the  curves  at  r/R  =  0.  40  have_scme  large  differences* 

In  particular,  for  240°  <  Y  <  360°,  the  ft,  ,  C,  ,  and  flt  Glauert 
coefficients  have  different  signs  for  the  two  advance  ratios. 

COMPARISON  OF  H-34  DATA  WITH  UH-1A  DATA 


The  outboard  sections  for  the  two  helicopter  blades  at  the 
lower  advance  ratios  have  certain  similarities  that  can  be  seen  by  com¬ 
paring  Figure  9  with  Figure  15.  However,  the  more  inboard  stations 
become  less  and  less  similar,  although  the  excursions  remain  roughly 
comparable.  At  the  higher  advance  ratios,  the  azimuthal  signatures 
appear  to  be  quite  different  except  for  the  very  low  frequency  content, 

COMPARISON  OF  TWO-DIMENSIONAL  AIRFOIL  DATA  WITH  ROTOR 
DATA 


It  is  clear  from  the  previous  discussion  of  the  results  that  the 
Glauert  coefficient  representation  of  pressure  differentials  is  a  sensitive 
indicator  of  the  shape  of  the  distribution.  Further,  strong  azimuthal  and 
radial  dependencies  are  indicated  in  Figures  7  through  18,  This  point 
is  further  reinforced  by  the  "hovering"  data  (  ,u.  =  0)  for  the  UH-1A 
shown  in  Table  X  where  the  following  ranges  are  noted: 
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Evidently,  the  division  of  circulation  between  the  C0  and  C,  distributions 
in  the  quantity  (C0t^C,)  varies  strongly  with  radius.  It  is  of  interest  to 
compare  these  data  with  the  two-dimensional  data  because  this  flight 
condition  should  have  the  smallest  unsteady  effects. 
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The  basis  for  the  comparison  was  equal  normal  force  coeffi¬ 
cients  and  equal  Mach  numbers.  Test  parameters  given  in  Reference  10 
were  used  to  generate  the  following  estimates  for  the  UH-1A  hovering 
condition: 


at  r/R  =  0.  40 


M  »  0.  30 


CL  *  1.12 


at  r/R  =  0.  85 


M  ~  0.55 

co 

CL  »  0.78 


at  r/R  =  0.  90 


0.  58 
CL  *  0.  72 


Section  data  curves  in  Reference  5  were  entered  with  these  values  and 
corresponding  geometric  angles  of  attack  were  found  to  be  12  degrees, 

7.  5  degrees, and  7  degrees  for  r/R  =  0.  40,  0.  85,and  0.  90,  respectively. 
These  values  for  o<  and  A^were  used  to  obtain  Glauert  coefficients 
by  interpolation  in  Table  VIII.  The  following  comparison  results: 
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The  agreement  between  the  Glauert  coefficients  for  the  two  cases  is 
surprisingly  good  in  view  of  the  various  estimates  required  to  make  this 
particular  comparison.  Of  course,  the  finite  span  effects  must  have 
influenced  the  results  at  r/R  of  0.  85  and  0.90.  Nevertheless,  the 
behavior  is  apparently  dominated  by  the  Mach  number  and  angle  of  attack 
(the  Reynolds  numbers  in  the  two  cases  for  the  two  sets  of  data  were  in 
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the  range  from  one  to  about  four  million  so  that  no  large  Reynolds 
number  effect  would  be  expected).  The  tentat;ve  conclusion  to  be  drawn 
here  is  that  for  rotor  stations  inboard  of  90  percent  radius,  most  of  the 
viscous  effects  are  similar  to  those  encountered  with  two-dimensional 
airfoils.  There  is  required,  however,  an  estimate  of  the  angle  of  attack 
for  such  usage  of  two-dimensional  wind-tunnel  data. 
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SERIES  REPRESENTATION  OF  THICKNESS  EFFECTS 


THEORETICAL  VALUES 

To  find  the  pressure  distribution  due  to  the  thickness  of  the 
NACA  0012  and  0015  airfoils,  it  is  necessary  to  cc/isider  the  effect  of  a 
round  nose  and  finite  trailing  edge  angles.  In  Reference  6,  it  is  shown 
that  the  pressure  due  to  the  thickness  of  the  airfoil  is  related  to  the 
surface  slope  of  the  airfoil  through  the  following  two  equations 


tx.  * 


sin 


#♦ 


e, 


-  t )  SLtt  9  +  £  Qn  sin.  nO  ^  ^ 


ft-Ml 


ft 


sin 


~JU 


iayl  T 


~  +  (18) 


where 


+  Z  8*  Cosn& 

=  the  thickness  function  of  the  airfoil 
t  = 

=  leading  edge  radius 
9t  =  trailing  edge  angle 
Cp  =  pressure  coefficient 
Bn  =  unknown  coefficient  to  be  computed 


For  NACA  four -digit  wing  sections,  the  thickness  U  is  given  by 
(Reference  7) 


y  *  5T  -  .lZbX  -  .351(0  X*  +  . 2 0+3 *3-.  7075%*] 


(19) 


where  X  ,  Tf  ,  and  T  are  nondimensionalized  by  the  chord  Zb,  and  X  - 
maximum  thickness.  The  leading  edge  radius  is  given  by 

^  =  I. !0/9  r* 

Differentiating  Equation  (19)  with  respect  to  %  ,  one  obtains 
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cL% 


«  sr 


f  .1484-5 

L  4% 


-  JZ I*  -  .7032  X  +  .  BSZ1  X  -.4-OUX 


(20) 


The  coefficient  of  5m  &  in  Equation  {17)  is  defined  to  be  8,  and  is 
given  by 


8. 


A 

Tf 


-t 


The  remaining  coefficients,  Bn  ,  in  Equation  (17)  are  determined  by 
standard  Fourier  expansion  techniques  which  depend  on  the  orthogonality 
of  the  sine  (and  cosine)  functions  to  give 


^  6  .  &t 

catTr  — 


sinZ  y  )  sm  n&  d.  9  ;  n&2  (21) 


Coefficient  values  presented  below  were  obtained  by  numerical 
integration  of  Equation  (21)  after  substitution  of  Equation  (20). 


NACA  0012 

HACA  0015 

B! 

-0.0686 

-0.1108 

b2 

-0.0304 

-0.0380 

B3 

-0.0090 

-0.0113 

B4 

-0.0070 

-0-0900 

B5 

-0.0070 

-0.0088 

B6 

-0.0047 

-0.0059 

«7 

-0.0040 

-0.0050 

B8 

-0.0023 

-0.0050 

B9 

-0.0011 

-0.0014 

B10 

0 

0 

COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  OF 
THICKNESS  EFFECTS 

Calculated  pressure  coefficients,  Cp  ,  due  to  thickness  are 
compared  with  experimental  results  at  zero  angle  of  attack  but  with 
different  Mach  numbers.  Figures  3  and  4  show  good  agreement  between 
computed  and  experimental  results. 
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PROFILE  DRAG  ESTIMATES 


THEORETICAL  ANALYSIS 


A  method  proposed  by  Helmbold  (Reference  11)  is  used  for 
estimating  profile  drag.  The  profile  drag  coefficient  for  each  surface  of 
the  airfoil  of  chord  length  c  has  the  following  form: 


(22) 


where 


CD 

% 


V, 

K, 

Ut 


profile  drag  coefficient 

Reynolds  number  based  on  chord  and  the 
free-stream  velocity 

position  of  boundary -layer  transition  from 
laminar  to  turbulent  flow 


local  tangential  velocity  at  the  outer  edge  of 
the  boundary  layer 


free-stream  velocity 

(&) 


U 2.5  R 


<u 


( 


3.1f 


V 


boundary-layer  momentum  thickness  at  the 
point  of  transition 


tangential  velocity  at  the  outer  edge  of  the 
boundary  layer  at  the  laminar-to-turbulent 
transition  point 


The  factor  K,  is  determined  by  the  condition  that  the  laminar  and  tur¬ 
bulent  momentum  thicknesses  are  equal  at  the  point  of  transition.  Hence, 
ftt  is  generally  not  the  same  on  the  upper  and  lower  surfaces. 

The  velocity  distribution  at  the  outer  edge  of  the  boundary 
layer  is  related  to  the  pressure  measurements  through  Bernoulli's 
equations,  i.  e. , 


-,«3 


(23) 
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or 


=  measured  pressure  coefficients. 

Equation  (23)  applies  to  both  upper  and  lower  surfaces  of  the 


References  4  and  5  provide  pressure  measurements  at  both 
upper  and  iower  surfaces,  but  only  pressure  differentials  are  given  in 
References  8,  9*  and  10.  Therefore,  the  pressure  coefficients  on  the 
upper  and  lower  surfaces,  Cp^  and  respectively,  are  constructed 
by  adding  the  experimental  thickness  effects  to  the  experimental  lifting 
effects.  Let 


where  Cp 
airfoil. 


where 

4>t  =  measured  pressure  coefficients  due 

to  thickness  (at  zero  angle  of  attack) 

=  velocity  due  to  thickness 

Then  again  using  Bernoulli's  equation  for  the  upper  and  lower  surfaces  of 
the  airfoil,  one  obtains 


-ft. 


(25) 


/•  _  tl  -  Ho  f  .  Wcx2 

»  f  p  v;  -  - 


(26) 


and 


=  H  ' 


(27) 
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where 


■fi 


pressure  at  upper  surface 
pressure  at  lower  surface 
velocity  due  to  camber 


For  a  given  set  of  and  £4  ,  the  velocity  due  to  thickness 
can  be  computed  by  Equation  (24),  while  Equation  (27)  provides  the  value 
of  v.c  .  Then  and  can  be  found  by  substituting  ut  and  tte  into 

Equations  (25)  and  (26),  respectively.  It  should  be  noted  that  near  the 
stagnation  point,  the  second  order  effect,  namely,  x  Ut/V, ,  has  to 

be  included.  Away  from  the  stagnation  point,  “*/ y ,  is  small  to  compare 
with  unity  and  Equation  (27)  reduces  to  a  simpler  form. 


The  position  of  the  transition  point  Xt  ,  on  the  upper  surface, 
is  assumed  to  be  located  at  the  point  of  minimum,  pressure.while  the 
transition  point  on  the  lower  surface  is  at  th#  leading  edge. 

The  boundary  layer  momentum  thickness  at  the  point  of  transi¬ 
tion  is  given  by  Reference  11  (Equation  37  in  Chapter  12): 


Equations  (22)  and  (28)  were  integrated  numerically  by  the 
trapezoidal  rule;  results  are  plotted  in  Figures  19  through  23. 


COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR 
TWO-DIMENSIONAL  CASES 

The  experimental  data  for  the  NACA  0012  airfoil  section 
presented  in  Reference  4  were  obtained  on  a  section  of  an  H-34  blade, 
and  the  drag  was  measured  in  two  independent  ways:  (1)  by  means  of  a 
force  balance  and  (2)  by  means  of  an  integration  of  a  wake  survey.  These 
results  are  indicated  by  the  curves  through  the  data  points  in  Figures  20, 
21,  22,  and  23.  One  of  the  interesting  aspects  here  is  that  the  two  methods 
of  measurement  give  quite  different  results.  The  percentage  errors 
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occurring  at  angles  of  attack  greater  than  4  degrees  are  exceptionally 
high.  In  general,  there  is  a  tendency  to  believe  that  the  wake  survey 
data  is  more  representative  of  the  two-dimensional  airfoil  than  is  the 
balance  data.  This  situation  occurs  because  the  wake  survey  can  be 
made  near  the  centerline  of  the  tunnel  and  can  be  relatively  characteristic 
of  the  two-dimensional  flow  case;  whereas,  the  balance  data  is  influenced 
by  side-wall  effects  and  by  the  usual  difficulties  of  developing  an  adequate 
drag  balance.  The  calculated  drag  as  a  function  of  angle  of  attack  and 
Mach  number  was  determined  from  the  pressure  distributions  and  was 
based  on  the  presumption  of  attached  flow.  Therefore,  the  calculated 
variation  cannot  exhibit  the  large  changes  that  occur  with  flow  separation 
at  angle  of  attack  near  the  stall  angle.  Nevertheless,  it  can  be  observed 
that  the  calculation  based  on  the  measured  pressure  distribution  agrees 
well  with  the  wake  survey  data  up  to  the  point  of  incipient  stall. 

The  corresponding  drag  coefficients  for  the  0015  airfoil 
calculated  from  Equation  (22)  are  compared  with  the  experimental  results 
from  Reference  5  in  Figure  24.  The  theory  overpredicts  the  value  of  the 
drag  at  the  lower  angles  of  attack  up  to  the  critical  Mach  number.  At  an 
angle  of  attack  of  8  degrees,  however,  the  measured  and  theoretical 
values  have  come  into  good  agreement  at  Mach  numbers  below  the  critical 
Mach  number.  The  trend  of  these  coetticients  with  angle  of  attack  was 
correctly  estimated.  Again,  the  theory  did  not  account  for  separation 
effects  nor  did  it  account  for  the  occurrence  of  local  shock  waves. 

ESTIMATES  OF  ROTOR  PROFILE  DRAG  DISTRIBUTIONS 

Since  the  drag  calculation  appeared  to  have  merit,  ar.  effort 
was  made  to  estimate  the  rotor  blade  drag  characteristics  on  the  basis 
of  the  measured  pressure  distribution.  These  results  are  summarized 
in  Figures  25  through  28.  In  Figure  25  there  is  presented  the  azimuthal 
variation  of  the  calculated  drag  coefficient  for  r/R  =  0.  40,  and  advance 
ratios  of  /x  =  0.  18  and  jul  =  0.  29  for  the  H-34  (that  is,  0012  ainoil 
section).  The  characteristics  of  the  curves  on  Figure  25  are  that  the 
drag  coefficients  remain  relatively  constant  for  0  <  Y  <  180  degrees. 

The  primary  difference  between  the  two  advance  ratios  was  the  large 
increase  experienced  in  the  vicinity  of  "Y  -  270  degrees  at  the  inboard 
section  (  r/R  =  0.40)  for  the  higher  jx  case.  On  Figure  26  there  are 
plotted  the  variations  of  drag  coefficients  at  the  two  outboard  stations, 
r/R  =0.  85  and  0.  90.  The  drag  coefficients  have  about  the  same  char¬ 
acteristics  as  those  noted  at  r/R  =  0.  40  except  that  the  amplitudes  of  the 
outboard  variations  are  somewhat  smaller. 

The  corresponding  drag  data  for  the  TJH-1A  are  shown  in 
Figures  27  and  28.  The  characteristics  previously  pointed  out  for  the 
H-34  data  are  also  in  evidence  in  the  UH-1A  data.  The  levels  for  the 
coefficients  are  very  similar  for  the  two  rotors,  and  the  variations  with 
advance  ratio  are  also  quite  similar. 
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CONCLUSIONS 


Both  flight  measured  rotor  and  two-dimensional  wind-tunnel 
differential  pressure  distributions  of  corresponding  airfoil  sections  were 
fitted  by  a  Glauert  series  representation,  and  the  individual  terms  of 
the  series  were  assessed  for  their  physical  significance.  Inspection  of 
these  data  and  corresponding  theoretical  results  lead  to  the  following 
conclusions: 


1*  The  theoretical  result  that  the  symmetric  airfoils 
require  only  one  Glauert  coefficient  to  represent 
the  pressure  differential  distribution  is  not  sup¬ 
ported  by  the  experimental  data  for  the  NACA 
0012  and  0015  since,  for  example,  appreciable 
camber  terms  were  found. 

i 

2.  The  Glauert -coefficients  representing  two- 
dimensional  experimental  data  experience  varia¬ 
tions  with  Mach  number  and  angle  of  attack  which 
are  not  readily  explainable.  These  variations  may 
involve  an  interaction  between  boundary  layer  and 
compressibility  phenomena. 

3.  Pressure  distributions  measured  on  rotors  in 
flight  as  well  as  measured  on  airfoils  in  wind 
tunnels  have  strong  similarities  but  also  some 
notable  differences.  Apparently,  the  viscous 
and  compressibility  effects  observed  on  two- 
dimensional  airfoils  are  largely  carried  over 

when  the  same  section  is  installed  on  the  rotor  blade. 

4.  The  characteristic  pressure  differential  components 
deduced  from  UH-1A  and  H-34  flight  data  are 
similar. 

5.  The  magnitude  of  the  higher  Glauert  components 
evident  in  the  rotor  pressure  distributions  suggest 
complex  pitching  moment  variations  and  concomitant 
control  loads. 

6.  The  pressure  distribution  analyses  indicated  that 
unsteady  effects  tend  to  be  concentrated  at  the 
inboard  stations  of  the  blade  on  the  retreating 
side  of  the  disk.  However,  no  attempt  was  made 
to  inspect  results  outboard  of  r/R  =  0.  90. 
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Rotor  section  drag  estimates  based  on  boundary 
layer  theory  appear  to  be  useful,  at  sections  where 
the  flow  does  not  separate.  Rotor  drag  distributions 
have  never  been  measured,  so  it  is  necessary  to 
rely  on  calculation  procedures  to  obtain  such 
estimates  analytically. 
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recommendations 


Based  on  the  results  of  this  study,  it  is  recommended  that: 

1.  Further  attempts  be  made  to  extract  additional 
information  from  the  available  helicopter  rotor 
blade  pressure  distribution  measurements. 

Finite  span  effects  and  Mach  number  effects 
are  among  these. 

2.  An  attempt  be  made  to  design  an  airfoil  which, 
when  operating  in  the  regime  of  interest  to 
helicopter  rotors,  will  introduce  negligibly  small 
control  loads.  Judicious  introduction  of  geometric 
camber  is  recommended.  In  particular,  the  camber 
components  of  the  Glauert  series  would  furnish 
required  characteristics.  It  is  believed  that  it  would 
be  necessary  to  account  for  viscous  effects  in  such  a 
design. 

3.  Attempts  be  made  to  relate  stalled  airfoil  data  in 
which  rate  effects  are  present  to  the  distributions 
measured  on  the  retreating  side  of  the  disk  on 
helicopter  rotors. 

4.  Transonic  airfoil  sections  be  designed  to  operate  in 
an  environment  deduced  from  the  available  rotor 
measurements  made  at  high  shaft  speed  and  high 
translational  speeds. 
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THEORETICAL  GLAUERT  COEFFICIENTS  FOR  NACA  0012  AIRFOIL  WITH 
TRAILING  EDGE  FLAP  (NEGATIVE  THREE  DEGREES  FLAP  ANGLE) 
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TABLE  TS. 


A: 


•7 


0.301 


0.40 


0.50 


0.60 


0.65 


0.70 


0.75 


oi 


A 


.0135 

.0586 

.0423 

.0127 

.0384 

.0432 

.0085 

.0408 

.0456 

.0089 

.0629 

.0647 

.0080 

.0801 

.0824 

.0054 

.0968 

.1031 

.0027 

.0948 

.0909 

.0001 

.0878 

.0192 


FIRST  THREE  GLAUERT  COEFFICIENTS,  AQ)  A,,  A2, 

OBTAINED  FROM  MEASURED  PRESSURE  DIFFERENTIALS 
ON  NACA  0032  AIRFOIL  SECTION 


4  6  8 


10 


12 


14 


16 


.0317 

.0863 

.0762 

.0301 

.0787 

.0875 

.0264 

.1059 

.1057 

.0157 

.1431 

.1462 

.0206 

.1339 

.1540 

.0116 

.1670 

.1743 

.00840 

.1372 

.1024 

.0039 

.1063 

.0315 


.0511 
.1008 
.1177 
.0508 
.0973 
.1338 
.0385 
.1542 
.1643 
.0279 
.1947 
.2101 
.0255 
.2193 
.2363 
.0172 
.1974 
.201 1 


.0741 

.1479 

.1438 

.0776 

.1360 

.1543 

Tom 

.2334 

.2567 

.0298 

.2313 

.2490 

.0756 

.2314 

.2014 

.0226 

.2224 

.1715 


.0980 

.1570 

.1731 

.loi  2 

.1455 

.1737 

.0490 

.2162 

.2412 

.0352 

.2274 

.1745 

ToiTo 

.2268 

.1344 

.0239 

.2327 

.1146 

.0168 

.2659 

.1002 


.1109 

.1641 

.2008 

T0739 

.2117 

.2201 

.0578 

.1923 

.1455 

'.0470 

.2056 

.0979 

.0401 

.2027 

.0696 

.0323 

.2120 

.0615 


.1167 

.1654 

.2128 

.0609 

,2135 

.1156 

.0372 

.2309 

.1138 

T040? 

.2102 

.0484 

.043? 

.2104 

.0579 

,042_5 

.1920 

.1733 


.0169 

.2618 

.0799 

.0388 

.2312 

.0726 

.0344 

.2395 

.0970 

.0324 

.2338 

.0582 

.0331 

.2881 

.1232 

7034? 

.2830 

.0774 


1 


0.8{ 


1 

‘2 


TABLE  I.  FIRST  THREE  NOND 1  MENS  1 ONAL  GLAUERT  COEFFICIENTS, 
A0,  Aj,  A2,  OBTAINED  FROM  MEASURED  PRESSURE 

DIFFERENTIALS  ON  NACA  0012  AIRFOIL  SECTION 


Mao  A 


.3)56 

.4116 

1.3688 

1.1768 

.9886 

1.0393 

.3983 

.4331 

1.2034 

1.1338 

1.3545 

1.2599 

.2927 

.3325 

1.4147 

1.3351 

1 . 5786 

1.3320 

.2200 

.1803 

1 . 5600 

1.6394 

1.6035 

1.6756 

.2355 

1.6687 

1.5291 

1.7157 

1.7582 

.1011 

.1228 

1.7978 

1.7544 

1.9152 

1.8310 

.0539 

.1091 

1.8923 

1.7818 

1.8137 

1.3299 

.0237 

.0679 

1.9953 

1.8642 

.0435 

.5521 

6 

8 

.5037 

.5005 

.9927 

.9991 

1.1590 

.9711 

.5108 

.5329 

.9784 

.9341 

1.3455 

1.0600 

.3333 

.1961 

1.3334 

1.6079 

1.4300 

.2228 

.2051 

1.5545 

1.5898 

1.6775 

1.7116 

.1886 

.1813 

1.6228 

1.6373 

1.7488 

1.4250 

.1482 

0.1689 

1.7035 

1.6623 

1.7357 

1.2817 

10 


.5552 

.8895 

.9805 


.5819 
.8363 
.  9984 


.3121 

1.3758 

1.5349 


.2415 

1.5169 

1.1642 


.2146 

1.5707 

.9307 


.8172 


12 


.5748 

.8504 

1.0406 


.4114 

1.1773 

1.2236 


.3751 

1.2499 

.9453 


.3166 

1.3668 

.0508 


.2835 

1.4331 

.4919 


.2337 

1.5327 

.4448 


.5851 

.8299 

1.0673 


.3632  .2513 

1.2736  1.4374 


.2168 
1.5664 
.3317  I  .3902 


.2916  .1867 

1.4168  1.6267 

.1 


.3898  .0956 


.3066  .1963 
1.3868  1.6074 
.1251  .4399 
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TABLE  31.  FIRST  THREE  GLAUERT  COEFFICIENTS,  A0,  A|,  A2,  OBTAIHED 
FROM  MEASURED  PRESSURE  DIFFERENTIALS  ON  NACA  0012  AIR¬ 
FOIL  SECTION  WITH  TRAILING  EDGE  FLAP  (THREE  DEGREES 

FLAP  ANGLE) 

Moo 

*1 

a 0 

0 

2 

4 

6 

8 

10 

12 

14 

16 

0.30 

*0 

-.0047 

.0030 

.0270 

.0445 

.0672 

.0884 

.1044 

.1217 

.0794 

*1 

-.0072 

.0131 

.0544 

.0926 

.1199 

.  1 429 

.1556 

.1535 

.1552 

*2 

-.0024 

.0408 

.0808 

.1237 

.1560 

.1888 

.2200 

.2226 

.0373 

0.40 

*0 

-.0054 

.0095 

.0268 

.0456 

.0693 

.0969 

.0859 

.0290 

.0446 

A1 

-.0084 

.0214 

.0667 

.1053 

.1283 

.1334 

.1806 

.2445 

.2119 

A2 

-.0071 

.0458 

.0922 

.1345 

.1651 

.1803 

.2280 

.1272 

.0644 

0.50 

*0 

-.0076 

.0088 

.0243 

.0387 

.0292 

.0452 

.0596 

.0371 

.0370 

A1 

-.0066 

.0292 

.0786 

.1380 

.2155 

.2203 

.1863 

.2210 

.2191 

A2 

.0093 

.0525 

.1069 

.1630 

.2572 

.2556 

.1445 

.1185 

.0864 

0.60 

*0 

~.0058 

I 

.0062 

.0136 

.0247 

.0238 

.0340 

.0445 

.0419 

.0353 

A1 

-.0163 

.0433 

.1131 

.1708 

.1976 

.2139 

.1972 

.1912 

.2334 

A2 

-.0008 

.0636 

.1383 

.2043 

.2196 

.1167 

.1175 

.0579 

.0925 

0.70 

*0 

-.0061 

.0019 

.0080 

.0141 

.0209 

.0222 

.0271 

.0347 

.0301 

A1 

-.0204 

.0617 

.1408 

.1844 

.2076 

.2193 

.2017 

.1976 

.280  4 

A2 

-.0021 

.0884 

.1727 

.2151 

.1772 

.1385 

.0721 

.0555 

.0934 

0.80 

*0 

-.0065 

-.0001 

.0045 

A1 

-.0460 

.0672 

.0932 

.0365 

.0172 

.0885 

J 
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TABLE  m.  FIRST  THREE  NONDIMENS  1 ONAL  6LAUERT  COEFFICIENTS,  A0,  Aj,  A2, 
OBTAINED  FROM  MEASURED  PRESSURE  DIFFERENTIALS  ON  NACA  OCI2 
AIRFOIL  SECTION  WITH  TRAILING  EDGE  FLAP  (THREE  DEGREES  FLAP 
ANGLE) 


*1 

a 

O 

*oo 

0 

2 

4 

6 

8 

10 

12 

14 

16 

0.3 

*0 

.5672 

.5512 

.4985 

.4900 

.5287 

.5530 

.5730 

.6133 

.5057 

A1 

.8657 

.8976 

1.0028 

1.0199 

.9427 

.8940 

.8541 

.7734 

.9886 

*2 

.2888 

2.8033 

1.4905 

1.3625 

1 . 2263 

1.1814 

1.2081 

1.1217 

.2375 

0.4 

*0 

.6036 

.4717 

.4456 

.4645 

.5195 

.5922 

.4876 

.1916 

.2965 

A» 

.7929 

1.0565 

1.1088 

1.0711 

.9611 

.8155 

1.0245 

1.6169 

1.4071 

h 

-.6750 

2.2617 

1.5314 

1 . 3687 

1 . 2368 

1.1023 

1.2938 

.8399 

.4277 

0.5 

*0 

.6965 

.3727 

.3819 

.3593 

.2)34 

.2910 

.3902 

.2513 

.2525 

*| 

.6069 

1 . 2547 

1 . 2362 

1.2813 

1.5733 

1.4180 

1.2196 

1.4974 

1.4950 

*2 

-.8500 

2.2522 

1.6810 

1.5133 

1.8778 

1.6452 

.9458 

.8028 

.5899 

I0 

.4153 

.2227 

.1935 

.2246 

.1677 

.2414 

.0311 

.3048 

.2323 

ft 

A  | 

1.1694 

1 . 5546 

1.6131 

1.5508 

1.6646 

1.5172 

1 . 3782 

1.3905 

1.5354 

1 

*2 

.0538 

2.2839 

1.9721 

1.8553 

1 . 4204 

1 . 3242 

.8212 

.4207 

0.7 

Bl 

.3759 

.0579 

J016 

.1329 

.  1940 

.1684 

.2115 

.2600 

.1769 

A! 

1.2483 

1.8841 

1.7968 

1.7341 

1.6120 

1.6632 

1 . 5770 

1.4801 

1.6461 

*2 

.1314 

2. 6987 

2.20  46 

2.0228 

1.7919 

1.6507 

.5638 

.4153 

.5485 

TABLE  um.  FIRST  THREE  GLAUFRT  COEFFICIENTS,  A0,  Aj ,  A2  AND  X0,  A,,  A2, 
_ OBTAINED  FROM  MEASURED  PRESSURE  DIFFERENTIALS  ON  NACA  0015  SECTION 


ES-’-.'i 


TABLE  IX.  FIRST  THREE  HOND IMENS I ONAL  GLAUERT  COE 
OBTAINED  FROM  H-34  FLIGHT  DATA  (0012  A 


TABLE  X.  FIRST  THREE  NOHD IMENS I OHAL  GLAUERT  COEFFICIENTS,  (T0,  C 
OBTAINED  FROM  UH-IA  FLIGHT  DATA  (0015  AIRFOIL  SECTION) 


Figure  3.  COMPARISON  OF  EXPERIMENTAL  PRESSURE  DISTRIBUTIONS 

(REFFRF.NCE  4)  WITH  THEORETICAL  RESULTS  AT  MACH  NUMBERS 
0.3  AND  0.4  TO  ILLUSTRATE  THICKNESS  EFFECTS. 


50 


Figure  «l.  COMPARISON  OF  EXPERIMENTAL  PRESSURE  DISTRIBUTIONS 

(REFERENCE  1)  WITH  THEORETICAL  RESULTS  AT  MACH  NUMBERS 
0.5  AND  0.6  TO  ILLUSTRATE  THICKNESS  EFFECTS. 
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x/c 


Figure  5.  SCHEMATIC  OF  UPPER  AND  LOWER  BOUNDS  FOR 
MEASURED  PRESSURE  DIFFERENTIALS. 
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54 


Figure  13.  COMPARISON  OF  THEORETICAL  NOND I  MENS  I OHAL  GLAUERT  COEFFICIENTS 
WITH  FLIGHT  DATA  FOR  UH-1A  AT  u  =  0.08,  r/R  =  0.40. 


60 


* 

Figure  17.  COMPARISON  OF  THEORETICAL  NOND I  HENS  I ON AL  GLAUEPT  COEFFICIENTS 
WITH  FLIGHT  DATA  FOR  UH-1A  AT  [i  =  0.26,  r/R  =  0.85. 
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— — EXPERIMENTAL  DATA,  AjVvf 
- EXPERIMENTAL  DATA,  C;  'n 


REGION  OF  LEADING  EDGE  SEPARATION  WITH  LAMINAR  REATTACHMENT 


=  angle  of  attack 

*■  =  LEADING  EDGE  RAD  I  US /CHORD  LENGTH 

R/  =  REYNOLDS  NUMBER  BASED  ON  LEADING  EDGE  RADIUS 


66 


a,  ANGLE  OF  ATTACK 


Figure  20.  COMPARISON  OF  MEASURED  AND  CALCULATED  DRAG  COEFFICIENTS 
FOR  MAC A  0012  AIRFOIL  SECTION  AT  MACK  NUMBER  0.3. 


SECTION  DRAG  COEFFICIENT 


Figure  25.  VARIATION  OF  CALCULATED  DRAG  COEFFICIENTS  WITH  AZIMUTHAL 
ANGLE  \p  FOR  H-34  AT  r/R  =  0.40,  fx=  0.18  AND  0.29. 


72 


Figure  27.  VARIATION  OF  CALCULATED  DRAG  COEFFICIENTS  WITH  AZMMUTHAL  ANGLE  ^ 
FOR  UH-1A  AT  r/R  =  0.40  AND  u  =  0.0,  0.08,  AND  0.26. 


APPENDIX  I 

LEAST-SQUARE  POLYNOMIAL  APPROXIMATION 


If  a  function  {  (X)  approximated  by  the  form 


{(X)  =  £  at  <pt  cx) 

r  0 


(29) 


is  to  hold  over  a  set  of  N+1  points  X0  ,  X,  , . ,  Xw  +  1  ,  where  N  in, 

=  unknown  constants,  and  fyCZ)  =  known  function  of  X  ,  and  the  sum 
weighted  square  error  is  to  be  a  minimum 


N 


Z 


-  Yru.nLYnu.rn. 


(30) 


where  -uX(Xi)  =  weighting  function  at  %i  ,  (  tJlXi)  =  I  for  all  J.  if  all 
the  data  are  of  equal  significance),  then  the  normal  equations  require  that 


or 


or 


3ar 

a / 


Z  md  [ f(Xz)  -£  a,  ■  *  0 , 

+  L  *.o  *  J 

L  <f>  (Xi)  f  { l%-)  -  Z  A.  *A7lA  «  0, 

L  *.0  *  *  J 

^  A/ 

L  i  CXi)  +  a,  l  -uj-(Xi)  <f>,(Zi  )  f 


(31) 


X  <f>r  (*;>  <f>nLXi)  =  t  Mi)  <f>,CXi)flX;) 

^*0  *'  -A 


aL-0 


(r  *  o,  i  , 
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Equation  (31)  is  a  set  of  n  +  1  linear  algebraic  equations  for  n  +  l 
unknowns  (i.  e. ,  a,  ,  a,  , . ,  ). 

Suppose  hS(Zj  =  1,  f(Z)~  kjp(&)  with  X  =  "  y  C05  ft  ,  A/+J  = 
total  number  of  pressure  measurements,  and 


X,  *  Re  cot  ■  +  Z,  sun.  k  ft 

k--t> 


Then  Equation  (31)  is  reduced  to 


f  ff  •£  t/tyfCO;)  ,  cr  =  o,  1 , 2,  ■  • ' ,  n) 

J  Vt*J 


or 


v  ft-  ft 

L  cot  f  Cot  -f 

JL 

7  SiM  cot  -f 

jL*o  *- 

N  -Q, 

X  sin  Z  Or.  cot 

X I'O  ^ 


L  SKaYi<9\  cot  -g  £  5 tn  n  ft  sin  ft.  •  • . .  £  sin.  n  ft  sm  nft 


y  %  ■  A 
Z  cot  —  5m  £► 


v*.*0 


7  sm  sin. 

<c*o 

7  sinZQ  stn<ft 


X  'fl 


V 


£  ft  • 

•  Z  eat  Sin.  n  ft 

T*  P  M 

JL*°  C“ 

// 

■  Z  Sin  ft  stn  n  ft- 
,  *  ^ 

A,  *0 

f/ 

7  stn.  ?.  ft.  stn  n  (ft. 

*-*  X  x 

x«o 

I 

I 


x  *0 


iZ. 


I 

2k 


A/  ft' 

Z  cot  -f  f«%) 

I  Stn  ft  ffftu 

«c  »© 

V 

£  stn  2ft  ff£.) 

X--0  ( 

I 
I 

1  J 

I  Z  sLnrL&.f(%) 
L  X  «  # 


(32) 
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It  should  be  noticed  that  all  the  pressure  differentials  (  W  +  I 
pressure  taps)  are  used  to  provide  n  +  1  Glauert  coefficients  (  A.  through 
fln  )  with  n  *  N  .  Equation  (32)  can  be  solved  by  the  standard  sub¬ 
routine  of  the  solution  of  the  *1  linear  simultaneous  algebraic  equations. 
Once  the  s  are  calculated,  quantities  like 
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